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We study quantum gravity with the Einstein-Hilbert action including the cos¬ 
mological constant on the Euclidean Einstein universe S^ x S^. We compute 
exactly the spectra and the heat kernels of the relevant operators on 5^ and use 
these results to compute the heat trace of the graviton and ghost operators and the 
exact one-loop effective action on 5' x 5^. We show that the system is unstable 
in the infrared limit due to the presence of the negative modes of the graviton and 
the ghost operators. We study the thermal properties of the model with the tem¬ 
perature T = (27rai)"^ determined by the radius a\ of the circle 5We show that 
the heat capacity Cy is well defined and behaves like ~ in the high temperature 
limit and has a singularity of the type ~ {T - Tc)~^, indicating a second-order 
phase transition, with the critical temperature Tc determined by the cosmological 
constant A and the radius a of the sphere S^. We also discuss some peculiar prop¬ 
erties of the model such as the negative heat capacity as well as possible physical 
applications. 
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1 Introduction 

The low-energy effective action in quantum field theory is a powerful tool that 
enables one to study the vacuum state of the theory. The low-energy effective 
action cannot be computed in the usual perturbation theory, and so to study it 
in the generic case, one needs new essentially non-perturbative methods. The 
development of such methods for the calculation of the heat kernel was initiated 
in our papers [1, 3] for a gauge theory in flat space, which were then applied to 
study the vacuum structure of the Yang-Mills theory in [4, 6]. These ideas were 
first extended to scalar fields on curved manifolds in [2, 5] and finally to arbitrary 
twisted spin-tensor fields in [8]. In [9] we applied these methods to study quantum 
gravity and Yang-Mills theory on any symmetric space. Further, we applied these 
methods to study the thermal Yang-Mills theory on product of spheres, such as 
5^ xS^ X 52 and5^ X 53 in [11, 12]. 

In the present paper we apply these methods to study the one-loop low-energy 
effective action in quantum Einstein general relativity in the Einstein Universe 
background at finite temperature. Erom the mathematical point of view, we com¬ 
pute the one-loop effective action for the Einsten-Hilbert action with cosmological 
constant on the background x S^. 

This paper is organized as follows. In Sec. 2 we introduce all the relevant 
operators for the calculation of the one-loop effective action in Einstein quantum 
gravity. We refer to the paper [9] for the details. In Sec. 3 we compute the heat 
trace coefficients. In Sec 4. we study the quantum gravity on 5^ and compute 
all relevant heat traces on 5^. We refer to the paper [12] for the details of the 
calculation of the heat traces on 5^ for any representation. In Sec. 5 we compute 
the heat traces and the effective action on 5 ^ x S^. Einally, in Sec. 6 we discuss 
the thermodynamic properties of the model. 


2 One-loop Einstein Gravity 

In this section we follow our previous work [9, 10]. The dynamics of the gravita¬ 
tional field parametrized by the Riemannian metric on a closed (compact without 
boundary) manifold (M,g) of dimension n is described by the ffilbert-Einstein 
action of general relativity, which (in Euclidean formulation) has the form 

f dxg^'H-^ + 2A) , 

M 


( 2 . 1 ) 
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where g = det^^y, G is the gravitational eonstant and A is the eosmologieal eon- 
stant. The elassieal vaeuum Einstein equations are determined by the first varia¬ 
tion of the aetion 

CO 1 

16;rGg-^^^-= R^'' - -Rg^^ + = 0 . (2.2) 

6g^lv 2 


In two dimensions the action is trivial 


S = {-Anx{M) + 2Avol(M)}, 
16;rG 


(2.3) 


where xi^) is the Euler characteristic of the manifold M and vol (M) is its vol¬ 
ume. Therefore, it does not have any extremal metrics; more precisely, in two 
dimensions every metric satisfies the Einstein equations with zero cosmological 
constant, 

1 

Rab - IjRgab, (2.4) 

and, therefore, the Einstein equations do not have any solutions for any K Q, 
which means that Einstein gravity in two dimensions is purely topological. 

Eor this reason, we restrict ourselves to n > 2. In this case the Riemann tensor 
can be decomposed as follows 


R^Kd = + - - - R5^‘'^c5'’\v 

cd n-1 ^ ^ (n - l)in - 2) ^ ^ 

where Cabed is the Weyl tensor. The norm of the Riemann tensor is then 

4 2 

n nabed _ /-* /^abed , r> nab r}2 

r<^abcdr^ - ^abed^ “I-“ 7-7T7- . 

n — 2 {n - l){n - 2) 

The solutions of the Einstein equations determine the Einstein spaces, 

2 

Rab — 1 

n-2 


and, therefore. 


2n 

R = -A. 

n-2 


In this case the Riemann tensor is 


nab 
« cd 


C‘^\d + 


(n - l)(n - 2) 




(2.5) 

(2.6) 

(2.7) 

( 2 . 8 ) 

(2.9) 
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with the norm 

D itabcd _ /^abcd , ^ A 2 /O 1 riA 

KabcdK - ^abcd^ + 7 -777-7777-^ • t^-tU; 

{n - l){n- ly 

The case of three dimensions is special. In this case the Weyl tensor is equal 
to zero identically, and, therefore, the Riemann tensor is fully determined by the 
Ricci tensor, 

( 2 . 11 ) 

Therefore, in particular, 

RabcdR‘^'’‘'‘' - ARabR'^” + R" = 0. (2.12) 

The Einstein equations take the form 

Rab = 2Ag,t (2.13) 

and, therefore, the curvature tensor of Einstein spaces is fully determined by the 
metric, 

(2.14) 

This means that the only Einstein spaces in three dimensions are the (locally) 
maximally symmetric spaces, the sphere 5^ for A > 0, the hyperbolic manifolds 
//-Vrfor A < 0, where E is a lattice in S 0^{\, 3); for A = 0 the only solutions are 
flat manifolds, like a torus . In any case, gravity in three dimensions is rigid, 
that is, it does not have any propagating degrees of freedom. 

Notice that the same invariant (2.12) plays a role in higher dimensions as well. 
In particular, in dimension n = A the integral of that invariant determines the Euler 
characteristic of the manifold 

XiM) = ^^dxg^l\RabcdR‘^>’^‘^-ARabR!^'’+R^) 

M 

^ sh I -^RabR^^’ + Ir"^. (2.15) 

M 

When the Einstein equations (2.7),(2.8), are satisfied the Ricci tensor is deter¬ 
mined by the metric. That is, for n = 4, 


Rab ^Sabt 


R = 4A, 


(2.16) 
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and the integral norm of the Riemann tensor is determined by the Euler eharacter- 
istic 

;^(M) = ^ J dx g^l^RabcdR‘^‘’^‘^ 

M 

= + (2.17) 

M 


Note that the Euler characteristic of Einstein spaces in four dimensions is positive 
definite. It is worth stressing that this disagrees with eq. (115) in [18]. 

The diffeomorphism invariance of the Einstein-Hilbert functional means that 
the metric carries some non-physical (gauge) degrees of freedom described by a 
vector field. In n dimensions a vector field has n independent components and a 
symmetric 2-tensor field has n{n -l- l)/2 independent components. Therefore, the 
gravitational field in n dimensions has 


n(n+l) n(n-3) 

N(n) = --- 2n = --- 


(2.18) 


degrees of freedom. This number is equal to N(4) = 2 in four dimensions as 
expected; however, it vanishes in three dimensions, N{3) = 0. In two dimensions 
it gives a meaningless result, N{2) = -1. We will compute the effective action 
in three dimensions below but one should realize that in three dimensions the 
Einstein gravity does not have any dynamics [19, 20]. 

One of the fundamental problems of quantum Einstein gravity is that the Eu¬ 
clidean Einstein-Hilbert action is unbounded from below, which leads to the diver¬ 
gence of the Euclidean path integral over all metrics. This divergence is concep¬ 
tual in nature and is much more serious than the usual ultraviolet divergence of the 
quantum field theory. It is well known [18] that under a conformal transformation 


gl,y = 


(2.19) 


where m is a smooth positive function on M, the action takes the form 
1 8 (n-l) 


16:7rG {n 


^ r 
-2) J 

M 




( 2 . 20 ) 


where Y is the Yamabe operator 


n — 2 
4{n - 1)' 


-R. 


y = -A + 


( 2 . 21 ) 
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The Yamabe operator is nothing but the conformally covariant scalar Laplacian. It 
is a self-adjoint elliptic partial differential operator with a positive leading symbol. 
The spectrum of such operator is real, discrete, and with finite multiplicities; it 
is bounded from below and unbounded from above. This shows that the action 
functional is unbounded from below. It is obvious that by keeping the metric g^y 
constant and taking the function co to be bounded and increasingly oscillating the 
action can be made arbitrarily large and negative. This is a well known conformal 
problem of quantum gravity. It has been suggested [18] that this problem can 
be avoided by deforming the contour of integration in the path integral over the 
conformal factor to make it purely imaginary, which will turn the action into a 
standard functional of quantum field theory. However, such an approach cannot 
be taken seriously. This is a major problem of Einstein quantum gravity and it 
remains open. A solution to this problem would require a modification of the 
Einstein-Hilbert action but we do not attempt to solve it in the present paper. 

The standard loop expansion of the Euclidean effective action has the form 

Y = S + ^Ed) + 0(n\ (2.22) 

where r(i) is the one-loop effective action. The one-loop effective action is de¬ 
termined by the graviton operator L 2 acting on symmetric two-tensor fields and 
the Eaddeev-Popov ghost operator Li acting on vector fields. In the Euclidean 
formulation the zeta-regularized one-loop effective action has the form 

r(i) = -^4«(0), (2.23) 

where 

^GR{s) = a^is)-2a,{s), (2.24) 

and and are the zeta functions of the operators Li and L 2 defined by 

00 

Uis) = d-^e-^^"GdO, (2.25) 

0 


where 

0l(O = Tr exp(-tL). (2.26) 

The renormalization parameter p is introduced to preserve dimensions and z is a 
sufficiently large infra-red regularization parameter, which should be set to zero 
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at the end of the ealeulation. Therefore, 

oo 

^gr{s) = ^ fdt r^e-^^"0GR(t), (2.27) 

^(s) J 

0 

where 

0GR(t) = 0L,(t)-20L,(t)\ (2.28) 

we will eall this invariant the heat trace of quantum gravity. 

The operators Li and L\ are determined by the second variation of the action 
and then by imposing some gauge condition on the metric fluctuation (see, for 
example, [15, 7]). The second variation of the action defines a second-order partial 
differential operator P acting on symmetric two-tensors by 

\enGg-^l \ Kp = , (2.29) 

dgf^vSga/} 2 


where 


pfjv,al3 _ _ _ g‘^l^g^^^^ A 

_gA;vv(Q'v^) _ 


_ ga(npv)P _ ^(iipv)a pRvpaPgRv 

+ (R - 2A). (2.30) 

In the minimal gauge the non-diagonal derivatives in both the graviton opera¬ 
tor and the ghost operator vanish and the operators take the form 


cd,ab 

2 


f ab 


^ga(Cgd)b _ ^gabgcd^ (-A + R- 2A) 

_2p(c\a\d)b _ ^a{cpd)b _ gb{cpd)a) pcd^ab ^cdpab 

-/*A - 


(2.31) 

(2.32) 


We should stress that the operator Li differs from the eq. (16.37) in [14]. 
The tensor 

1 


■pcd,ab 


ga(Cgd)b _ -g^^bgCd 


(2.33) 
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here is the metrie in the spaee of symmetrie tensors. It is easy to see that it is 
positive definite in the subspaee of traeeless symmetrie tensors but it is negative 
definite in the eonformal (sealar) seetor. This is exaetly the problem of the eonfor- 
mal mode in quantum gravity diseussed above. Following the standard approaeh 
[18, 13, 14, 15] we simply assume that it ean be fixed somehow by some physieal 
arguments and proeeed as follows. We faetor out this metrie from the operator 
L 2 to define the graviton operator L 2 and the ghost operator Lj in the eanonieal 
Laplaee-type form 



Lj - A + Qj , 

(2.34) 

where the potentials for both operators are [9] 


(20% = 

Tya 

-K b , 

(2.35) 

{Q2r\b = 

-2R\a\) - 26^\aR^\) + R^'^gab + -^g^'^Rab 

n-2 

-T-^^S^'SabR + 6\a6\){R - 2A). 

{n - 2) 

(2.36) 


We should stress here that the endomorphism Q 2 does not eoineide with the eq. 
(16.78) in [14]. 

It is well known that the heat traee of Laplaee type operators has the asymp- 
totie expansion as t ^ 0 

00 

(2.37) 

k=0 

where Bk{E) are the so-ealled Hadamard-Minakshisundaram-DeWitt-Seeley coef¬ 
ficients (or simply heat trace coefficients) of the operator L. This means that the 
function 0(t) has similar asymptotic expansion as t ^ 0 

00 

@GR{t) - {Antr'^Yj^^Cu, (2.38) 

k=Q 

where 

Cu = Bu{L2)-2Bu{LQ. (2.39) 

It is easy to find the dependence of the effective action on the renormalization 
parameter; by integrating the equation 

d 

= -^gr(O) , 


( 2 . 40 ) 
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we get 

r(i)Ou) = r(i)Ouo) - log j ^gr{0) . (2.41) 

This enables one to study the high-energy asymptotics of the effective action as 

jX ^ CO. 

For the Laplace type operators the value of the zeta function at 5 = 0 is deter¬ 
mined by a specihc heat trace coefficient 


Therefore, 


^l(O) = 


(4;7r) for even n, 

0 , for odd n. 


4’gr(0) - 


(47r) '‘^^Cn/ 2 , for even n, 
0 , for odd n. 


in particular, in four dimensions, n = 4, 


^CRiO) = (47Tr^C2. 


(2.42) 


(2.43) 


(2.44) 


3 Heat Trace Coefficients 

We will need the heat trace coefficients Bq, Bi and B 2 for the operators Lj and Li. 
They have the following well-known form [17, 7] (we neglected the inessential 
total derivatives here which do not contribute to the global invariants) 


Bo(L) 

II 


(3.1) 


M 

(1 



Bi(L) 

II 

-e). 

(3.2) 


M 

\ ^ 



B2(L) 

II 

fl 2 

\-Q 

2^ 

-\RQ + ^nabn‘^'’ 

0 IZ 



M 


+/ ( J_7?2 + _Lp pabcd _ J_^ I (3 3) 

\72 180 180 // 
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Here / is the identity endomorphism and "Rah is the eurvature of the spin eonnee- 
tion of a tensor field realizing a representation of the spin group defined by 

Kv = (3.4) 

where Ijab are the generators of the spin group Spin(n) satisfying the eommutation 
relations 

~ ~§ac^bd "I" Sbc^ad Sad^bc ~ Sbd^ac- (3.5) 

For the veetor representation the identity and the generators have the form 



Hifd 

= 5^d, 


(3.6) 


^(\),abf d 

= '^5‘^[agb]d, 


(3.7) 

Therefore, tr /j 

= n and 





5o(Li) : 

= nvol (M). 


(3.8) 

Also, we have 

trQi 

= -R, 


(3.9) 

and, therefore. 

BdU) = f 

dx g^'^\:{n + 6)R. 

0 


(3.10) 


kJ 

M 




Further, we eompute 








(3.11) 

and 

^^R{l),abR{l)^^ 

D ryabcd 

- -^abcd^ 


(3.12) 

We also have 

trceo' 

II 


(3.13) 

and, therefore. 





52(Li) = J 

M 

D uabcd ,90 

dX g ■! ^gQ RabcdR + ^gQ RabR + 

n + 12 1 
72 J 

(3.14) 
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Next, we eompute the heat traee eoefheients for the operator L 2 . For the tensor 
representation the identity and the generators have the form 


ihr^cd = (3.15) 

{^i2),abf'cd = ^5^\agbMc5^\). (3.16) 

First, we have tr /2 = n{n + l)/2, hence, 

5 o(L2 ) = in(n+l)vol(M). (3.17) 

Now, we introduce the following endomorphisms 


(Fi)^",, 

_ -pc d 
~ ^ {a b)") 

(3.18) 

(V2r‘'ab 

= 

(3.19) 

{v^r^b 

= R^^'gab, 

(3.20) 


II 

(3.21) 

iv^r^b 

II 

(3.22) 


Then the endomorphism Q 2 has the form 

Q2 = Q2 + {R-2K)l2, (3.23) 

where 

Q 2 = -2yi - 2y2 + V 3 + (3.24) 

n -2 {n — 2 ) 

Now, by using the traces 


tryi = --R, 

(3.25) 

try2 = ^(n+1)^. 

(3.26) 

trya = R, 

(3.27) 

tr V 4 = R, 

(3.28) 

tr ys = n. 

(3.29) 

tr Q 2 = -nR, 

(3.30) 


we compute 
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trViVi 

= 

3 

D nabcd 

-^Kabcdt<- , 

(3.34) 

trV2V2 

= 

^-{n + 2)RatR‘^^ + ]^R\ 

(3.35) 

trVaVa 

= 

R\ 

(3.36) 

trV4V4 

= 

R\ 

(3.37) 

trVsVs 

= 

n\ 

(3.38) 

trViy2 

= 

^ D T?ab 

--i<abK , 

(3.39) 

tryiy3 

= 

D nab 
l<ab^ , 

(3.40) 

trViV4 

= 

Ty nab 

t<abt< , 

(3.41) 

IrViVs 

= 

R, 

(3.42) 

try2V3 

= 

n nab 
^abJ< > 

(3.43) 

tvV2V4 

= 

n nab 

Kabts , 

(3.44) 

try2V5 

= 

R, 

(3.45) 

try3V4 

= 

nRabR^\ 

(3.46) 

try3V5 

= 

nR, 

(3.47) 

try4V5 

= 

nR, 

(3.48) 


we obtain 

tx{Q2f = 3RabcdR‘'’’^‘' + ~ (3.49) 

n—2 n—2 
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By using these results we get 




8n + 4 


n-2 

-2n(n - l)RA + 2n(n + 1)A^. 


P r>ai> 1 
Kabi< + 

2 


n^ - 5n^ + Sn + 4 

2(n - 2) 


R^ 


(3.50) 


We introduce yet another endomorphism 


(TpcY'^ab = S^\pg,yaS‘^b). 

By using 

tr Tp,r^ = -l(n + 2)S^\p6^\^ 

we get 

tri:^2iab^t.2r = -4(n + 2)d[^[,d''],], 

and, therefore, 

trRaiabR(2Y‘’ = -{n + 2)RabcdR^'’^Y 

By using these results we get 

B2(L2) = { dx g''^[a,RabcdR"^^'' + a2RabR‘‘'' + a,R^ + n^A + 72 ^^ 

= 29.-480), 

-n^ + 181^2- 1438n +720 
360(n - 2) ’ 

25n^ - I45n^ + 262n + 144 
144(n - 2) ’ 

= ^n(7 - 5n), 
o 

= n{n + 1). 


where 


J 

M 


a\ 


a2 


as 


71 

72 


(3.51) 

(3.52) 

(3.53) 

(3.54) 

(3.55) 

(3.56) 

(3.57) 

(3.58) 

(3.59) 

(3.60) 
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Finally, by using these results we eompute the coefficients Ck, (2.39), 


Co 

1 

2 

n{n - 3)vol (M), 

(3.61) 

Cl 

= J 

dx g^'^ < 

|-Y^(5n2 - 3n + 24)7? + n{n + 1)a| , 

(3.62) 


M 



C2 

- J 

^dxg^l^\ 

PlRabcdR^^^‘' +P2RabR‘^^ +PiR^ + + 72 A") 

9 


M 



(3.63) 

where 


Pi 

= T4:(n'-33n + 540), 

JoU 

(3.64) 



Pi 

-n^ + 185n2 - 1806n + 1440 

360(n - 2) 

(3.65) 



P3 

25n^ - 149^2 + 222n + 240 

144(n - 2) 

(3.66) 


These results are also different from the ones given by eqs. (16.79)-(16.81) of 
[14]. 

We will need these coefficients for n = 4. By using the above results we obtain 
in this case 


Bo(^i) = 4vol(M), 

Bi(Li) = J dxg^'^^R, 

M 

= r dx g^l^[-^RabcdR"''"‘^ + ^RabR"^ + 

J ® \ 180 90 9 j 


(3.67) 

(3.68) 


B2(LQ 


(3.69) 
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5o(L2) = 

lOvol(M), 



5i(L2) = 

f dxg^l^< 

[-Y^ + 20 a| 

5 


M 

V ) 


B2{L2) = 

f 

19 

D rycibcd 

—Kabcd^ - 

55 


(3.70) 

(3.71) 


36^ “ T' 


M 


Therefore, we obtain the total coefficients (2.39) 
Co = 2vol(M), 

23, 


r 

( 


Cl = ( dxg^<U-—R + 20A\, 


(3.72) 

(3.73) 

(3.74) 


M 


C2 


f dx [ ^Rabc 4 R‘^'’‘^‘^ - ^RabR‘“’ + ^R^ - ^RA + 20A^ 


M 

1696 


45 


90 


36 


n^XiM) + fdx g^'^ [-^RabR^^ + ^R^ - ^RA + 20A2^ 


90 


36 


M 


(3.75) 


Finally, when the Einstein equations are satisfied these coefficients in four dimen¬ 
sions take the form 


Co = 2vol(M), 

(3.76) 

Cl = -^Avol(M), 

(3.77) 

Cl = n^xiM) ^A\o1 (M). 

(3.78) 

This gives the quantity 


106 29 A2 


^gr(0)= vol(M), 

45 40 

(3.79) 


determining the scaling properties of the model. It is worth stressing that this 
result coincides with the eq. (4.23) of [13] and disagrees with the results of [18], 
eq. (79), and [16], where the second coefficient is 657/540 instead of -29/40; in 
particular, it is not positive definite contrary to the claim of [18]. 
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4 Heat Traces on 


4.1 Reduction to Irreducible Representations 

In this section we compute the effective action and the relevant heat traces on the 
3-sphere of radius a. We define the dimensionless cosmological constant by 

A = a^A. (4.1) 


The curvature in the orthonormal frame has the form 

j^ah 




0 




cd ~ 

R = 


Here and below Sabc is the three-dimensional Levi-Civita symbol. 
The volume of the sphere is 


(4.2) 

(4.3) 

(4.4) 


vol(5^) = 2n^a^ 

and the Euclidean classical Einstein-Hilbert action on S ^ is equal to 

S = ^a(-3 -t a^A ). 


(4.5) 


(4.6) 


Note that the classical action is bounded from below and attains a minimum equal 
to 


2 GVa 

at the radius determined by the cosmological constant 

ao = 


(4.7) 


(4.8) 


so that classically the dimensionless cosmological constant is equal to 1, To = 1 • 
We will compute the heat trace of the Laplacian for the unit sphere 5 ^ by setting 
a = 1; the trivial dimensional factor a can be easily restored at the end of the 
calculation by replacing t tja^. 
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Let 

{WKd = ^g‘^‘’gcd (4.9) 

be the projection to the scalar representation and 

n 2 = / 2 -no (4.10) 


be the projection onto the space of traceless symmetric tensors. Note that the 
projection 112 acts as identity in the subspace of traceless symmetric tensors. In 
three dimensions, the dimensions of these subspaces are 



II 

II 

p\ 

(4.11) 


trflo = 1, trn 2 = 5. 

(4.12) 

This is consistent with the dimension of the general irreducible representation 
labeled by an integer j. 

Then the potential terms are 

tr A = 27 + 1 . 

(4.13) 

Qi 

= - 2 / 1 , 

(4.14) 

Qi 

= (4 - 2d)n2 - (2 + 2T)no . 

(4.15) 


We will reduce the calculation of the heat traces of the operators Li and L 2 to 
the calculation of the heat trace on the unit sphere S ^ of pure Laplacians Aj acting 
on irreducible representations j, 

&j{t) = Tr exp(tAy). (4.16) 

First of all, we immediately see that since the endomorphism Qi is constant, we 
have 

exp(-tgi) = e^‘Iu (4.17) 

and, therefore, the heat trace of the ghost operator is 

0L,(O = e2^0i(O . (4.18) 

We also have a similar formula for the operator L 2 , 


&L 2 (t) = Tr exp(-t22) exp(tA). 


(4.19) 
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However, the general tensor representation eontains the irredueible represen¬ 
tation with j = 2 (traeeless symmetrie two-tensors) and the sealar representation 
with J = 0 (traee). The spaee of symmetrie tensors deeomposes eanonieally into 
the direet sum of the traeeless tensors and sealars with the eorresponding projee- 
tions II 2 and Hq. It is easy to see that 

expC-tQz) = e^-™*Il2 + (4.20) 

Therefore, the heat traee of the graviton operator takes the form 

^ _ ^4 21) 

4.2 Heat Trace for Irreducible Representations 

Beeause the graviton operator neatly splits, we only need to eompute the heat 
traees for Laplaeians in irredueible representation j for integer j. This heat traee 
ean be computed by using the heat kernel diagonal for the Laplacian on the 
unit sphere given by the eqs. (6.16) of our paper [12]. To get the heat trace we 
have to multiply the heat kernel diagonal by the volume of the sphere S ^ equal to 
vol (5^) = and by the dimension of the representation j equal to ( 27 - 1 -1). This 
gives 


©.(;) = :^^-3/2^.[A;+ihi] ^ ^expj- 




n n 


n=-oo \n\<j 

Following [12] we introduce the function 


-p^t\ 1 - Ip^t - 




(4.22) 


Q(0 = ^ 


exp - 


2 2 

n n 


(4.23) 


which can be expressed in terms of the Jacobi theta function 

f2(O = 03(O,e-"'/^), 

and satisfies the following duality relation 



-03(0,0- 

n 


(4.24) 


Q(0 = 


(4.25) 
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This enables us to express the function 0y in terms of the function Q. as follows 


0y(O = _^r^/2gW+i)+i] ^ _ 2p^t)Q.{t) - ltQ:{t)]. (4.26) 


\r\<j 


By using the duality relation, the function 0y takes the form 


0/0 


,t[;U+l)+l] 


00 ^ 


-tin^+fp') 


(4.27) 


n=-oo 1^1 <j 

Finally by using the obvious equation 


,/ j 




0 , 


EE("'-''T' 

fi=-j n=-j 

we get the heat trace of pure Laplacian in the irreducible representation j 

00 r 

0/0 = Zi “ ^1) 

n=j+\ '- 


(4.28) 


+ ^ 2{n^ - /// exp [-t [n^ + - j(j + D - l]) ■ (4.29) 

In particular, the eq. (4.29) gives the eigenvalues and their multiplicities of the 
pure Laplacian acting on an irreducible representation J of S U(2). It is labeled by 
two integers n and p such that 


0 < p < j < n. (4.30) 

The eigenvalues are given by 

= n^ +p^ - j(j + 1) - 1, (4.31) 

and their multiplicities are 

dn,oi-Aj) = n\ (4.32) 

for /i = 0 and 

dn,,{-Aj) = 2{n^-p\ (4.33) 

for 1 < /i < j. 
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The minimal eigenvalue of the Laplacian -Ay is 

T„,i„(-Ay) = j (4.34) 

with multiplicity Jmin(-Ay) = (j + 1)^. In particular, this means that all Laplacians 
-Ay for 7 > 1 are positive and the scalar Laplacian Aq is non-negative, it has the 
obvious constant zero mode. 

We will need the functions ©o, ©i and ©2, 

0o(O = ^r^^^e^[Q(t)- 2 tQ.'(t)] 

00 

= Yj (4.35) 

n=i 

©i(0 = ^ 2(1 - 2t)e^^] Qit) - 2t [e^^ + 2e^^] Q'(0) 

00 

= Yj + 2 - l) , (4.36) 

n=2 

©2(0 = + 2(1 - 2t)e^‘ + 2(1 - 80^^'] f2(0 

-2t[e'^‘+ 2e^‘+ 2e^*] n'(t) 

00 

= Yj + 2 (n^ - 1) e-f(”"-6) + 2(„2 _ 4)e-f(”"-3) j. (4.37) 

n=3 

It is worth noting that the contribution of /i = 0 and // = 1 for 7 = 1 corresponds 
to the decomposition of the vector fields 

(Pf,=Aj; + Vf^cr, (4.38) 

where is the transversal (divergence free) vector, and the contribution of 7/ = 0, 
/i = 1 and /i = 2 for 7 = 2 corresponds to the decomposition of the trace-free 
symmetric tensor fields 

- ^g^yAcr, (4.39) 

where ip-^^ is the transversal (divergence free) tracefree tensor and cr is a scalar. 
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4.3 Heat Trace of Quantum Gravity 

We introduce the trace-free tensor part and the scalar part of the graviton operator 

L2 


Lf = noL2no, (4.40) 

if = Ii2L2li2. (4.41) 

Now, by using the eqs. (4.18) and (4.21) we compute the eigenvalues of the 
operators L\ and L 2 

= /l„,^(-Ai) -2 = n^ + 1 / -5, n>2, yu = 0,1, (4.42) 

d„,^(Lf) = d„,o(-Ao) - 2 - 2d = _ 3 _ 2d, n > I, (4.43) 

dn,n(E^2^) = d„_^(-A2) -I- 4 - 2d = -I-- 3 - 2d, n > 3, yu = 0,1,2. 

(4.44) 

The minimal eigenvalues are 

d^i„(Li) = -1, (4.45) 

d^in(Lf) = -2-2d, (4.46) 

= 6-2d. (4.47) 


Notice that the minimal eigenvalue of the ghost operator Ly is always negative, 
the minimal eigenvalue of the conformal sector of the graviton operator is 
negative for d > -1 and the minimal eigenvalue of the graviton operator in the 
traceless tensor part is negative for d > 3. That is, the graviton operator is 
positive only for negative cosmological constant when d < -1. 

Next, by using eq. (2.28) and (4.21) we get 

0 ^^(?) = e^-^^^^'^^@2{t) + e^™'®o{t)-2e^^®i{t). (4.48) 

We can write this either in terms of the function Q 

= [e^'^‘ {2e^‘ + 2e^' + 2^“') - 2e^' - 4e^'] Q'{t) (4.49) 

+ [e^^‘ (2e^‘ + 2(1 - 2t)e^‘ + 2(1 - 80^“') - 2e^' - 4(1 - 2t)e^'] Q(o|, 


0gr(O 
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which is useful in the ultraviolet limit as t ^ 0, or in the speetral form 
+ 4g(-l+2i)( _ _ ^2 

e^^* + 2{n^- l) ^-'(”'-2) + 2{n^ - 

-2n2e-'(«"-5) - 4 - i) (4.50) 

which is useful in the infrared limit as t ^ oo. 

When the classical Einstein equations are satisfied, that is, when d = 1, the 
heat trace simplifies to 

0 Gi?(O = _ 2t)e‘^^ + 2(1 - %t)e‘] Q(t) - 2t [-2e^* + 2/j Q'(t) 

OO 

= -4e‘-12 + J^ (-2 (n^ - l) + 2(n^ - . 

n=3 

(4.51) 

The infrared properties are described by the limit t ^ oo. By using the spectral 
representation of the heat trace we immediately get 

ScRit ) = - % e ‘ - 12 + 6>(e(22-6)^^_ (■4 52) 

The exponential growth of the heat trace indicates the presenee of the negative 
modes. 

It is instructive to study the asymptotics of the heat trace as t ^ 0. By using 
the asymptoties of the funetion as t ^ 0 

Q.{t) ~ 1, Q'(0 ~ 0, (4.53) 

we obtain 

0gr(O = {Co + tCi + t^C2 + 0(t^)| , (4.54) 

Co = 0, 

Cl = 2;r2(12d - 30), 

C 2 = 2;r2(12d^ - 24d - 3). 



where 


(4.55) 

(4.56) 

(4.57) 
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These coefficients coincide with the coefficients Ck given by the general formulas 
(3.61)-(3.63) in three dimensions, n = 3. Notice the absence of the constant term 
here. This is the feature of three-dimensional quantum Einstein gravity since it 
does not have any dynamics, that is, the number of degrees of freedom is equal to 
zero. 


5 Heat Traces on xS^ 

5.1 Reduction of Heat Traces 

In this section we study Einstein quantum gravity in the physical four-dimensional 
Einstein Universe. Since we would like to study the thermal effects at the same 
time, we consider the four-dimensional Riemannian manifold M = x with 
a circle S ^ of radius a\ and a sphere 5^ of radius a. So, all indices in this section 
are four-dimensional, that is, they run over 1,2, 3,4. 

Let h^b be the projection tensor on S ^ and q‘^b be the projection to S ^ so that 

S\ = q\ + h\, (5.1) 

and 

q\q\ = q\, h\q\ = 0, (5.2) 

h\ = 3, q\ = \. (5.3) 

Also, we introduce the Levi-Civita tensor Sabc on S ^ such that 

Sahcd^d = 0, (5.4) 

and 




(5.5) 

P pdec 

^abc^ 


(5.6) 

^ ^dbc 
^abc^ 

= 2h\. 

(5.7) 

Then the curvature is 

II 

2 

(5.8) 

II 

jhab-> 

a^ 

(5.9) 

R = 

6 

a^ ' 

(5.10) 
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The volume of the manifold M = 5' ^ x 5 ^ is 

vol(M) = An^aia^, (5.11) 

and the elassieal action is equal to 

2 

S = —aia(-3 + a^A). (5.12) 

2G 

Thus, the potential terms (2.35), (2.36), of the operators Li and L 2 are 

a\Qirb = -2h\, (5.13) 

aHQ2r‘’cd = {6-2A)6\c6\)-2h^\,h^\)-h^%d 

(5.14) 


We need to compute the heat traces of the Laplace type operators Lj = -A+Qj 
(2.34) on M = 5' x S^. We note that since the potential terms are constant we 
have 

exp(-tL,) = exp(-t2,)exp(tA), (5.15) 

and also 

exp(tA^‘^^') = exp(tA^‘)exp(tA^'). (5.16) 

Therefore, the heat traces ©l/O of the operators Lj can he computed as follows 


&Ljit) = 0 






where 

e^\t)= 

with Q(0 defined by (4.23), is the heat trace on the unit circle S ^ and 


Ql](t) = Tr exp(-tLf) 


(5.17) 


(5.18) 


(5.19) 


is the heat trace on the unit sphere S^. 

The heat trace 0^ (0 on 5^ was computed in our paper [12]. We consider 
a tensor representation of spin j of the spin group Spin(4) with generators 
satisfying the algebra (3.5). Recall that Spin(4) = SU{2)x S U{2). Therefore, the 
matrices 


1 
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satisfy the algebra (no summation over j!) 

[G(j)i,G(j)k\ = -s‘ikG(j)i (5.21) 

and form a reducible representation of the group SU(2); with the Casimir operator 
(no summation over j!) 

^ 0 ') = (5.22) 

We also dehne the matrix 

Gu)iy) = (5.23) 

where y = (y') is a unit vector. Let / be a real-valued function of x = (x') 6 R^. 
Let x' = ry‘, where r = |x| = -\/x'x; and y = (y') is the unit vector such that |y| = 1. 
Of course, the unit vector y lies on the unit sphere 5^ in R^. We introduce the 
average over the unit sphere S ^ of functions in R^ by 

(/) (^) = ^ J" dys2firy); (5.24) 

the integration goes over the unit sphere S ^ with the appropriate induced metric 
on S^. 

Then the heat trace of the Laplace type operator Lj = -A + Qj on the unit 
sphere has the form [ 12 ] 

0^(0 = ^ exp [-t(Glj^ + Qj - /y)] Sj (t), (5.25) 

where S j{t) 


Siit) 


E 


exp - 


2 2 

n n 


oo 

J ^ 


2r^-2 


2 2 

n n 


(exp[2rVtG(j)(y)]) . 

(5.26) 


5.2 Generators 

For the vector and the symmetric 2-tensor representation (3.7) and (3.16) the gen¬ 
erators have the form 


(GmiTd = sf,, 
(G^DiT^d = 2ey^cy\), 


(5.27) 

(5.28) 
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so that 



(G(i)(y))% = Si\y\ 

{G(2){y)r\d = 

(5.29) 

(5.30) 

We compute the Casimir operators 



= -2h^d, 

= -eh^fch^^d) + Ih^^Kd - Aq^fchf^d) 

(5.31) 

(5.32) 

and the sums 



(Gfi) + Qifd 

1 

II 

(5.33) 

iGl) + Qitcd 

= (6 - 2X)6\c5\) - + h‘^'’h,d 

-h^^cd - q‘^%d - - 2q^\cd. 

(5.34) 


5.3 Algebra of Constant Symmetric Endomorphisms 

First of all, for the vector representation we immediately obtain 

Lemma 1 

exp + Qi)] = li-H + e‘^‘H , (5.35) 

where H is the matrix of the projection H = Off)- 

To compute this exponential for the tensor representation we need to do some 
algebra. We dehne the following basis of endomorphisms acting on symmetric 
two-tensors in four dimensions 


rah 
^ cd 


(5.36) 

A^‘’cd 

= h^%ch^\). 

(5.37) 

nah 

cd 

II 

(5.38) 

/^ah 

C cd 


(5.39) 

jyih 
!-> cd 

= q‘'''hcd. 

(5.40) 

rpab 
^ cd 

= q^\ch^\). 

(5.41) 

rpab 
r cd 

II 

(5.42) 

First, we note the identity 

A + 2E + F = I, 

(5.43) 
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so that 

E = ^{l-A-F). (5.44) 

Of course, I is the identity. We compute the squares of these endomorphisms 


= 

= 

= 

= 

E^ = 
E^ = 

and their products 

AB = B, 
AC = C, 
AD = 0, 
AE = 0, 
AE = 0, 
BC = 3C, 
BD = 0, 
BE = 0, 
BE = 0, 
CD = B, 
CE = 0, 
CE = C, 
DE = 0, 
DE = 0, 
EE = 0, 


A, 

(5.45) 

35, 

(5.46) 

0, 

(5.47) 

0, 

1 

2^’ 

(5.48) 

(5.49) 

F, 

(5.50) 


BA = B, 

(5.51) 

CA = 0, 

(5.52) 

DA = D, 

(5.53) 

EA = 0, 

(5.54) 

EA = 0, 

(5.55) 

CB = 0, 

(5.56) 

DB = 3D, 

(5.57) 

EB = 0, 

(5.58) 

EB = 0, 

(5.59) 

DC = 3E, 

(5.60) 

EC = 0, 

(5.61) 

EC = 0, 

(5.62) 

ED = 0, 

(5.63) 

ED = D 

(5.64) 

EE = 0. 

(5.65) 
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Next, we define the following endomorphisms 


Pi = 

A 

(5.66) 

Pi = 

3 

(5.67) 

P 3 = 

2E, 

(5.68) 

Pa = 

F, 

(5.69) 

T = 

C + D, 

\{Pa-P2-T), 

(5.70) 

X = 

(5.71) 

n± = 

i (P 2 + P 4 + X). 

(5.72) 


By using the algebra of these endomorphisms one ean prove 

Lemma 2 The endomorphisms Pi , P 2 , P 3 and P/^form a set of orthogonal projec¬ 
tions satisfying 

P- = Pi, (5.73) 

PiPj = 0, if i ^ f (5.74) 

and 

Pi+P 2 + ^ 3+^^4 =/. (5.75) 

The dimensions of the corresponding subspaces are determined by the traces 


trPi = 5, trP 2 = 1, trP 3 = 3, trP 4 = 1. (5.76) 


Of eourse, the total dimension of the spaee of symmetrie two-tensors in four 
dimensions is 


5 + 1 + 3 + 1 = 10. (5.77) 

Lemma 3 1. The endomorphism X satisfies the equations 

XPi = PiX = P 3 X = XP 3 = 0 . (5.78) 

(P 2 + P 4 )X = X{P 2 + Pa) = X. (5.79) 

= P 2 + P 4 . (5.80) 

trX = 0, (5.81) 


It has the eigenvalue 0 with multiplicity 8 and simple eigenvalues - 1 , + 1 . 
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2. The endomorphisms n+ are the projections to the eigenspaces of X corre¬ 
sponding to the eigenvalues +1. They satisfy the equations 



ni 

= n±, 

(5.82) 


n_n+ = 

n+n_ = 0 , 

(5.83) 



= ±n±, 

(5.84) 


n+(P2 + ^4) = n+. 

n_(P 2 + Pa) = n_. 

(5.85) 


trfl 

± = 1, 

(5.86) 

Proof The projections P, act on the matrices C and D by 



PiC = 0 , 

CP, = 0, 

(5.87) 


PiC = c. 

CP2 = 0, 

(5.88) 


P3C = 0 , 

CP 3 = 0, 

(5.89) 


P4C = 0 , 

CP 4 = C, 

(5.90) 


PiD = 0, 

DP, = 0, 

(5.91) 


PiD = 0, 

DP2 = D, 

(5.92) 


P3D = 0, 

DP3 = 0, 

(5.93) 


P4D = D, 

DP4 = 0. 

(5.94) 

and, therefore, 


PiT = 0, 

TP, = 0, 

(5.95) 


PiT = C, 

TP2 = D, 

(5.96) 


P3T = 0, 

TP3 = 0, 

(5.97) 


PaT = D, 

TP4 = c. 

(5.98) 

so that 


P2T + TP2 = r. 

P4T + TP4 = T. 

(5.99) 

Also, we have 


T^ = XP2 + P 4 ). 

(5.100) 


By using these equations one can prove all the equations of the lemma. 
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Since the matrix X is orthogonal to the projections Pi and P 3 , it has an obvious 
eigenvalue equal to zero with multiplicity 8 = 5 + 3 making it essentially two- 
dimensional. It acts nontrivially only on subspaces spanned by projections P 2 and 
P4, which are both one-dimensional. Sinee it is obviously traeeless, the sum of its 
eigenvalues is equal to zero. It is easy to see that it has two non-zero eigenvalues 
±1. This follows from the eqs. (5.80). 

The matrices n+ are the eigenprojections corresponding to the eigenvalues ±1; 
this follows from the eqs. (5.84). 

We prove the following 

Lemma 4 

exp j-t[Gf 2 ) + Qi]] = [Pie^' + Pae"^' + n_ + U+e~^‘]. (5.101) 

Proof. We have 

G^ 2 ) + 22 = (2 - 2/1) 7 - 4Pi + 2X, (5.102) 

Therefore, 

exp {-t[G^ 2 ) + 22 ]) = exp {-t (2 - 2A)} exp (4tPi) exp(-2tZ) (5.103) 
We eompute 

exp(4tPi) = Pie^' + P 2 + P 3 + P 4 . (5.104) 

The only thing left to compute is the exponential exp(-2tX). By using 

= P 2 + P 4 , = Z, (5.105) 

we get 

exp(-2tZ) = Pi+ Ps + U_e^‘ + U+e-^‘. (5.106) 

This finally gives the eq. (5.101). 

5.4 Algebra of Symmetric Endomorphisms on 5 ^ 

Let y = (y') be a unit vector orthogonal to q^b, that is, satisfying y“ = We 

introduee two matriees 


b — y b ■> 


( 5 . 107 ) 
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and 

lb — b y yb • 

The square of the matrix Z is equal to 

(5.108) 

Z^ = -P, 

(5.109) 

and the matrix P is obviously a projection so that 


P^ = P, 

(5.110) 

PZ = ZP = Z, 

PH = HP = P, 

(5.111) 

(5.112) 

and 

trZ = 0, trP = 2, trH = 3. 

(5.113) 


First, we prove 


Lemma 5 The exponential of the matrix G(i)(y) is 

exp[2rG(i)(y)] = f- P + cos(2r)F + sin(2r)Z (5.114) 

with the trace 

tr exp(2rG(i)(y)) = 2 + 2cos(2r). (5.115) 

Proof. This follows from the faet that 


G(i)(y) = Z (5.116) 

and the eq. (5.109). 

Next, we introduee the following endomorphisms anting on symmetrie two- 
tensors 


'db 

cd 

N 

II 

(5.117) 

db 

' cd 

— Zj ((-Zj 

(5.118) 

rdb 

cd 

— 7^'^ tA’) 

~ ^ (c^ d)t 

(5.119) 

rdb 

cd 

II 

(5.120) 

db 

cd 

II 

(5.121) 

db 

cd 

= P^^’Pcd, 

(5.122) 

db 

cd 

= P^\cd, 

(5.123) 

'db 

cd 

= g‘^'’Pcd. 

(5.124) 
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We compute the traces 


liK = 0, (5.125) 

trL = -1, (5.126) 

trM = 5, (5.127) 

trA^ = 3, (5.128) 

trW = 0. (5.129) 

We need to compute the algebra of these endomorphisms. First, we have 

+ (5.130) 

= M, (5.131) 

= N, (5.132) 

= N, (5.133) 

LM = L, (5.134) 

ML = L, (5.135) 

KM = ^{K+W\ (5.136) 

KL = -W, (5.137) 

NM = N, (5.138) 

NL = L, (5.139) 

KN = VK (5.140) 

KW = ^(-N + L). (5.141) 


We prove the following 

Lemma 6 The exponential of the endomorphism G( 2 )(y) has the form 

exp[2rG(2)Cy)] = y(r)I + p(r)M + v(r)N + A{r)L + p{r)W + x{r)K, (5.142) 
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where 


y(r) = 

1, 

(5.143) 

yu(r) = 

2 cos(2r) - 2, 

(5.144) 

v(r) = 

1 3 

- cos(4r) - 2cos(2r) + 

(5.145) 

A{r) = 

i i cos(4r). 

(5.146) 

dir) = 

sin(4r) - 2 sin(2r), 

(5.147) 

x{r) = 

2 sin(2r). 

(5.148) 

with the trace 



tr exp[2rG(2)(y)] = 4 + 4eos(2r) + 2cos(4r). 

(5.149) 

Proof. We note that 

G(2)(y) = 2K. 

(5.150) 

Let 

J{r) = exp{4rK). 

(5.151) 

It satisfies the differential equation 



drJ = 4KJ 

(5.152) 

with initial condition 

m = I. 

(5.153) 

We decompose it according to 



J = yl + pM + vN + dL + r]W + xK. 

(5.154) 


Then by using the algebra of the matriees M, N, L, W, K we have 
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Therefore, the coefficients of this expansion must satisfy the differential equations 


drj = 0, (5.156) 

drP = -2x, (5.157) 

drv = -27], (5.158) 

drA = 27] + 2 x, (5.159) 

drT] = 2/i + 4v - 4/1, (5.160) 

drX = 2/i + 4y, (5.161) 

with the initial conditions 

y(0) = 1, yu(0) = v(0) = 4(0) = ?7(0) = x{U) = 0. (5.162) 


The solution of this system gives the result (5.143)-(5.148). Now the trace can be 
easily computed. 


5.5 Group Averages 

Next, we need to compute the group averages (5.24) of the functions given by 
(5.114) and (5.142). Thus, we need to compute the averages of the polynomials. 
We prove 

Lemma 7 The averages of the monomials are 



1 , 

0 , 


1 

2 m 


^(' 1*2 . . . ghk-lhk) 


(5.163) 

(5.164) 

(5.165) 


Proof. The first two equations are obvious. To prove the eq. (5.165) we 
consider the Gaussian integral 


r 


dx 


• ■ ■ X 


^2n 


K 


312 


• • • ^'2n-U2ii) 

n! 2 ^” 


(5.166) 


By changing the variables here by x’ = ry' and using the integral 

oo 


r 


dr r^"^^e 


-4 


r- (2n + 2)! 

vTT- 

(n+l)!22«+3 


(5.167) 


we get eq. (5.165). 
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Corollary 1 Let h be the projection onto the three-dimensional subspace V = 
of R^ and y = (yf be a four-dimensional unit vector lying in V. Then the eq. 
(5.165) is modified as follows 

//I . ..y^2k\ ^ \ (5.168) 

Ik + I 


We define the characters of an irreducible representation j of 517(2) by 

aM) = tr; {exp[2rG(y)(y)]) , (5.169) 


where tr j is the trace in the irreducible representation j. For an irreducible repre¬ 
sentation j the average of a group element over the 5^ is given by eq. (5.55) of 
[ 12 ]: 

(exp[2rG(y)(y)]) = ^ cosi2pr)Uj, (5.170) 

^ \iA<i 


so that 

xM) = E cos(2/ir). (5.171) 

\r\<j 

in particular, 


Xi(r) = 1 -I- 2cos(2r), (5.172) 

Xiir) = 1 -I- 2 cos(2r) -I- 2 cos(4r). (5.173) 


Using the averages of the monomials calculated above we obtain 


(Z) 

= 0, 

(5.174) 

{P) 

= 

(5.175) 

(M) 

= .ip,). 

(5.176) 

(N) 

7 2 

15^^ ^3^^’ 

(5.177) 

(L) 

= 1(2^2-7^0, 

(5.178) 

(1^) 

= 0, 

(5.179) 

(K) 

= 0. 

(5.180) 
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This allows us to calculate the group averages of the exponentials 
<exp[2rG(i)Cy)]) = h-^H + ^//eos(2r), 

<exp[ 2 rG( 2 )(j')]) = ^ [1 + 2 eos(2r) + 2 cos(4r)] Pi + P 2 + Pa 

[1 + 2eos(2r)] P 3 . 


5.6 Heat Trace of Operator L\ on 5' ^ 

To compute the functions S j we will need the integrals 


r 


dre ^ 


eos 


{ip Vtr) 


—00 

00 


r 


dre-''" 


cos 


(ip ^^tr) , 




By using these integrals we obtain 


f 


dr e ^ \ 2r 


In^n^ 


cos( 2 rVt) = V^l 1 - 2 t 


2n^n^ 


and we hnally obtain from (5.26) 


Suit) = Yj exp(- - 


2 2 

71 n 


1 - 


2n^n^ 


Ii--H\ + \\-2t- 


2n^n^ 


(5.181) 

(5.182) 

(5.183) 

(5.184) 

(5.185) 



(5.186) 
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The heat trace is calculated now hy using (5.25) and (5.35) 

j — oo 

<U) = ^r’'v|;exp| 

n=-oo 

2n^n^ 


xtr 1 - 


2 2 

71 n 


+ 1 - 2t 


t 

t 


n=-oo 

2n^n 


h-H + -He 




2 2 

71 n 


(5.187) 


X 1- 


{e^ + e^') + 2\\-2t- 


27T^n^ 


At I 


The heat trace can be written in terms of the function Q, 

ef^ (t) = [e^‘ + e‘ + 2{l- 2t)e^‘] Q(t) - 2t [e^^ + e^ + 2e^^] Q'(o|. 

(5.188) 

This heat trace has the asymptotic expansion as t ^ 0 


0 iy^) = (4 + lOt + 13t^ + 0(t^)). 

This can be put in the spectral form by using the identities 

2 /- 3/2 

Q(t) - 2tQ'(t) = —— y n^e-‘"\ 


s/n 

2t3/2 


neZ 


{l-2t)Q{t)-2tQ.'{t) = ^Yin^ -l)e 

yn ^ 


(5.189) 

(5.190) 

(5.191) 


neZ 


We get 




The asymptotic behavior of the heat trace in the limit t ^ oo is 

elY) = 4e‘ + 1 + Oie-^‘). 


(5.193) 
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5.7 Heat Trace of the Operator L 2 on 5 ^ 

Now, by using the integrals (5.183) and (5.184), we eompute first 


dr e '' 2r^ - 


1 - 


(5.194) 


dr g-'' - 


eos(4rVt) = V^(l-8t- ^^^ "^^(5.195) 


and then by using (5.26) we find 




(5.196) 




e“' + -Pi l-8t- 


Further, by using (5.101) and the algebra of symmetrie endomorphisms we get 


exp[-t(Gp) + exp 


(5.197) 


(ri / 27:'^n^\ 2 , 2 / 

X) - 1- \e^‘ + - \l-%t 

I [5 \ t j 5 \ 

1 / 27d-n^ \ -2t 2 / 

"3 '-— r "3 '-'^ 


+ - 1 - 2t 


Pi 


— e-^^ P3 + (n_ + 11 

t '' u 


Finally, by taking the traee we obtain from eq. (5.25) the heat traee 




X 1 


+ 2 1 - 2t 


+2 1 -2t 


e' + (3 - 16t 


+ 1 


e-^‘)\. (5.198) 
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We can rewrite this in terms of the function Q. as 

= -^r^/V'^'|[e3? + 2(l-20e^' + e' + (3-160e“' + 2(l-20e“2' 

+e-^‘]Q.{t) - 2t [e^* + le^‘ + e' + 3e-‘ + le-^‘ + Q'(o|. (5.199) 

To second order in t, the exponentials may be expanded 

[lO + (-26 + 20A)t + (35 - 52A + 20A^)t^ + 0(t^)]. (5.200) 
By using the identities (5.190)-(5.191) and 

2/3/2 

(3 - 160f^(t) - 6tf2'(0 = ^ y (3n^ - 8)e-^”\ (5.201) 

id 

the heat trace can be rewritten in the spectral form 

0^y) = y e-'”"{n^e^^ + 2 {n^ - l)e^‘ + n\' + {3n^ - 8)^“' 

+2(n^ - l)e~^‘ + 

The asymptotic behavior of the heat trace in the limit t ^ oo is 

0^J(/) = e2-i^[i + 0(g-3?)j (5.202) 


6 Effective Action 


Now, by using (5.17), (5.18), the heat trace of quantum gravity (2.28) on 5 ^ x 5^ 
takes the form 




( 6 . 1 ) 


It will be convenient to separate the asymptotic behavior at t ^ 0 


_ ^ ^ naia" _ / t \ t , 

0Gi?(O = T——1^( —1> 


4 f 


a^ 


(6.2) 
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where 

W{t) = (0 - 20^; (O). (6.3) 

Note that the funetion W depends also on the radius a through the dimensionless 
eosmologieal eonstant X = a^A. The asymptoties of the funetion IT as t ^ 0 are 

lT(t) = Cq + Cit + C2t^ + 0{t^\ (6.4) 

where the eoeffieients ct are eomputed from eqs. (5.189) and (5.200) 

Co = 2, (6.5) 

Cl = -46 + 20/1, (6.6) 

C2 = 9-52d + 20dl (6.7) 


The eoeffieients differ from the eoeffieients Q, (2.39), by the volume factor 
vol (5^x5^) = 4;7r^aia^ and a uniform factor a^^. The asymptotics of the function 
IT as t ^ oo are 

W{t) = \e^^^ - 8c' - 14 + 0{e-^‘) + . (6.8) 

sfn ' ' 


Now, following [12] 


the one-loop effective action can be presented in the form 

r(i) = -f-(;eiog^ + o}, (6.9) 

S a [ pI J 


where 


O 


Rgr (0 


00 



c '^0 2 + ^Ci + 2//oj t + (c2 + Ci/^o + dt) ’ 


2 2 4 4 

C 2 - z a c\ + z a , 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 


z is an infrared regularization parameter, and po is an arbitrary renormalization 
parameter. 

The total effective action including the classical term in the one-loop approxi¬ 
mation is 

T = ^aia(-3 + d) - |>Slog^ + o| + 0{fi\ (6.13) 

2G % a { pI ] 


We neglect the terms of order fi^ and set ^ = 1. 
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7 Thermodynamics 


The effective action is a function of two variables, T = Y{a\,a), where a\ is the 
radius of the circle S ^ and a is the radius of the 3-sphere 5^. The temperature T is 
determined by the radius of the circle ai by T = 1 KlnaQ and the spatial volume 
V of the system is the volume of the sphere S^, equal to V = 2n^a^. We introduce 
a dimensionless temperature 


ai 


so that the the temperature is 

T = :^. (7.2) 

2 na 


Then for a canonical statistical ensemble with fixed T and V the free energy F is 
determined by the effective action T by 


X 

F = TT = — r. 
2 na 


(7.3) 


By using the results of the previous section for the effective action we obtain the 
free energy 

^ ^ ^ ^ 

This enables one to compute all other thermodynamic parameters of the gravi¬ 
ton gas such as the entropy 



dF 

S =-= -2nadxF, 

dT 

(7.5) 

the energy 


E = F + TS = F - xdxF, 

(7.6) 

the pressure 

dF 1 dF 

dV 6n^a^ da ’ 

(7.7) 


and the heat capacity at constant volume 

dE d^F 

Cv = — = -T-;— = -2naxdlF. (7.8) 
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We see that the classical term and the renormalization term in the free energy 
(7.4) do not depend on the temperature; therefore, the entropy and the heat capac¬ 
ity do not depend on those terms. Therefore, the entropy and the heat capacity at 
constant volume are given by the derivatives of the function O, 

5 = - 5 ^ 0 , 

8 

c, = Jjra;* 

By changing the integration variable 1 a^t we rewrite the function O as 

oo 

O = J je-'^"‘^"[a[xh)wit)-RGRit)]. (7.11) 

0 

Differentiating the function O with respect to x, we get 

oo 

= 2xJ je-^‘‘'^'Q'(x^t)W(t), 

0 
oo 

lQ'(xh) + 2x^tQ"(xh)jW(t) 

0 

We will need the asymptotics of the function Q(t) obtained in [11]. We have 
as t ^ 0 

Q(t) = l + 2exp|-yj + 0(e-^^'/^) , (7.14) 

and as t ^ oo, 

Q(t) = -^\t^'^+ lt^'^e-'+ 0(e-^*]\ . (7.15) 

Because of the asymptotic behavior of D' as t ^ 0, the integrals converge at 
t ^ 0. The function W{t) increases exponentially at infinity with exponent e^ or 
e^^\ the function O has a singularity for < max|^,2A|. We may then view 
the maximum of these parameters as analogous to Aqcd, an infrared cutoff below 
which our analysis ceases to describe this system. 


(7.12) 

(7.13) 


(7.9) 

(7.10) 
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When z is taken to zero, the integrals (7.11) and (7.13) do not eonverge, and 
so we do not examine the free energy or entropy past this point. However, as we 
will see later, beeause the asymptotie behavior of Q! (t) + 2 tQ.'' (t) as t ^ oo is 
proportional to e~‘, the heat eapaeity may eonverge even in the limit z ^ 0. We 
deeompose W aeeording to 

W(t) = (e^^‘ - Se‘ - 14) + V(t) , (7.16) 

yn ' ' 


where the funetion V is exponentially small as t ^ oo. We may then split the 
integral for 5^0 into four parts: 

dl<!> = h+l 2 + h+l4, (7.17) 

where 


h 

h 

h 

h 


0 

00 

0 

00 

0 

00 

2 r (x^t) + 2xhQ" (^t)]V{t). 


(7.18) 

(7.19) 

(7.20) 

(7.21) 


0 


Notiee that the funetion in the integrand is exponentially small at infinity, 
namely, 

n'(t) + 2tQ'\t) = -^(2t - 3)e-' + 0(e-^') (7.22) 

yn 

Therefore, the integrals I 3 and I 4 eonverge for any x. The integral U eonverges 
only for x^ > 2A - a^z^ and the integral I 2 converges for > 1 - a^z^. Allowing 
the infrared cutoff z to go to zero, the integral I\ converges only for x^ > 2A and 
the integral h converges for > 1. 

Therefore, all of the integrals converge at high temperature but the heat capac¬ 
ity has a singularity either at the temperature Xc = V2T (for positive A > \ 12) ox 



1. G. Avramidi and S. J. Collopy: Quantum Gravity in Einstein Universe 


43 


2 ii Xc = 1 if 4 < 1/2 (including the case of negative cosmologieal constant A). 
Recalling that A = a^A and x = 2naT, this defines the critical temperature 


below which the system will undergo a phase transition. Notice that the smallest 
value of the eritical temperature is 



(7.24) 


The phase diagram of the graviton gas for the positive eosmological constant has 
the form illustrated on the graph Figure 1. 

We also study the high temperature limit as x ^ oo. The asymptotic behavior 
of the combination of derivatives of (7.22) implies that the high temperature 
limit corresponds to the limit of t ^ oo. We hnd the limit of U through Ij, replacing 
the Q functions by their leading asymptotics and integrating: 


h 


h 


h 


_16 

TTX' 

128 

nx ’ 
224 

nx 


(7.25) 

(7.26) 

(7.27) 


The integral /4 is evaluated by changing variables t t/x^ and using y(0) = 
VF(0) = 2; we get 

h'-Ax^v, (7.28) 

where v is the eonstant dehned by the integral 


V = 


/ I, 

-{q: {t) + 2tQ:'{t)]. 

0 


(7.29) 


The integral h dominates in the high temperature limit and determines the heat 
capacity per volume 


y 


4na^ 


= 2n^vT^ . 


(7.30) 
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Figure 1: Phase diagram of the graviton gas 


The dependence is characteristic of the photon gas and, as has been found in 
our previous paper [12], of the gluon gas as well. 

Next we study the behavior of the heat capacity near the critical temperature. 
By using eq. (7.22) and setting z = 0 we get 


/2 


32x^ 

n 



dt te 


-t{jr-2A) 


256x^ 
n 



dt te-‘^’^-^'>. 


We obtain, as x ^ 


h 


32A 



(7.31) 

(7.32) 


n' 


(7.33) 
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and as X ^ r 


1 9R 

h - - {x-\r\ (7.34) 

n 

The critical exponent of (-1) is indicative of a second-order phase transition. 

The temperature at which the phase transition occurs depends on the value 
of the cosmological constant. In the case that A > (2a^)“\ the phase transition 
occurs at the temperature Tc = sjAKln^) and if A < {2a^)~^ (also if A < 0) the 
phase transition occurs at the temperature Tc = \/{2na). The asymptotics of the 
heat capacity near the critical temperature are: if A > (2a^)~^ then as T Tc 

03/2 

C,- a^A^^^ (T - TcT^ , (7.35) 

n 

with Tc = ■\jA/(2n^). Since the heat capacity at constant volume as a function 
of temperature grows at infinity, this means that the heat capacity must have a 
minimum at some temperature T] > Tc. Further, if A < {2a^)~^ then as T Tc 

C, ~ (T - Tc )-', (7.36) 

na 

with Tc = \l{2nay, in the case that A = (2fl^)“' we have 

c, (r-r,)-i, (7.37) 

na 

with Tc = l/(27ra). This means that the heat capacity must vanish at some tem¬ 
perature T 2 > Tc (see Figure 2.). 


8 Discussion 

It is well-known that the gravitational action is unbounded from below and unsta¬ 
ble. The primary goal of this paper was to study the quantum gravitational field 
restricted to a set of manifolds which have an action that is classically bounded 
from below, and then examine how one-loop quantum effects disturb that stability. 
In order to calculate the one-loop effective action exactly, it is necessary to study 
a spacetime with a great degree of symmetry. We studied the thermal Einstein 
universe S ^ x S ^ with non-zero cosmological constant, varying the model only 
with respect to the radii, ai and a, of the circle S ^ and the sphere respectively. 
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Figure 2: Heat capacity as a function of temperature 


This spacetime is off-shell, so, strictly speaking, it does depend on the gauge of 
the quantum field. However, we used the generally accepted minimal covariant 
De Witt’s gauge in which all operators become Laplace type. 

We computed the exact trace of the heat kernels of all relevant operators, which 
enabled us to calculate the one-loop effective action exactly. The lowest value of 
the of the graviton operator can be chosen to be positive by adjusting the cos¬ 
mological constant, but the ghost operator always yields a negative eigenvalue, 
indicating an unstable mode for any radius of the Einstein universe. This may 
indicate a problem with the gauge condition and requires a detailed further study. 

We also studied the thermal properties of the model. We found that while the 
free energy and entropy are ultraviolet divergent, the heat capacity is well-defined 
even in the infrared limit. In the high-temperature limit, the heat capacity of the 
graviton gas has a dependence which is typical of a photon gas, and has also 
been found in our previous paper [12] to be consistent with a gluon gas. 

We also computed the asymptotics of the heat capacity near the critical tem¬ 
perature and found that the heat capacity has a branching singularity ~ (T - TcT^ 
at a finite critical temperature Tc given by (7.23). 

In the case of negative or small positive cosmological constant, A < (2a^)“\ 
the system exhibits a rather anomalous peculiar behavior with the negative heat 
capacity due to the presence of the unstable mode of the ghost operator. It is com- 
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mon in bound gravitational systems to have negative heat eapaeity. For instanee, 
the temperature of a black hole decreases as heat is added to it. The fact that the 
heat capacity changes sign at some temperature T 2 indicates that the system has a 
minimum internal energy at that temperature. 

It is interesting to play with the minimal value of the critical temperature given 
by (7.24). If we substitute the observed value of the cosmological constant, A ~ 
10“^^m“^, then the minimum critical temperature is approximately 

= ( 8 . 1 ) 

One can speculate that if the universe cools below the critical temperature, it is 
likely that some degrees of freedom would be frozen leaving a cosmic background 
thermal graviton radiation with temperature 

The techniques used in this paper are very general. It would be interesting 
to extend this model to higher-derivative quantum gravity or to supergravity, in 
which the one-loop action vanishes on-shell. 
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